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The Kolmogorov-Zakharov spectrum predicted by the Weak Turbulence Theory remains elusive 
for wave turbulence of flexural waves at the surface of an thin elastic plate. We report a direct 
measurement of the nonlinear timescale Tjvl related to energy transfer between waves. This time 
scale is extracted from the space-time measurement of the deformation of the plate by studying the 
temporal dynamics of wavelet coefficients of the turbulent field. The central hypothesis of the theory 
is the time scale separation between dissipative time scale, nonlinear time scale and the period of the 
wave (Td » Tnl >> T). We observe that this scale separation is valid in our system. The discrete 
modes due to the finite size effects are responsible for the disagreement between observations and 
theory. A crossover from continuous weak turbulence and discrete turbulence is observed when the 
nonlinear time scale is of the same order of magnitude as the frequency separation of the discrete 
modes. The Kolmogorov-Zakharov energy cascade is then strongly altered and is frozen before 
reaching the dissipative regime expected in the theory. 



PACS numbers: 46.40.-f,62.30.+d,05.45.-a 



I. INTRODUCTION 

The nonlinear coupling among a large number of waves 
causes an energy transfer between them and yields a state 
of wave turbulence. This phenomenon plays an impor- 
tant role in the atmosphere/ocean coupling through sur- 
face waves forced by the wind, in energy transport and 
dissipation by internal waves across the ocean, energy 
exchanges by Rossby waves in the atmosphere, energy 
dissipation by Kelvin waves in superfiuid turbulence or 
confinement in tokamak plasmas among many examples. 
The weak turbulence theory (WTT) aims at developing 
a theoretical framework for such turbulence of dispersive 
waves [J-Q- It relies on two major assumptions: weak 
nonlinearity and an asymptotically large size of the sys- 
tem. Because of the weak nonlinearity, energy exchanges 
between waves are slow and restricted to resonant waves 
that conserve energy and impulsion. The wave amplitude 
is slowly modulated by the coupling and an energy cas- 
cade operates from the large scales of the energy input 
to the small scales at which it is dissipated. The theory 
relies fundamentally on the following scale separation 



T < T NL < T t 



(1) 



where T = 2tt /lj is the period of the wave, Tnl the char- 
acteristic time of the nonlinear exchange of energy in the 
cascade and Td is the characteristic time of dissipation. 
Both latter times usually depend on u>. Thanks to this 
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scale separation, an asymptotic calculation can be per- 
formed to predict the statistical properties of wave tur- 
bulence [l|-[3( . It predicts in particular an energy cascade 
in scales sharing many similarities with the Kolmogorov 
cascade of hydrodynamical turbulence (thus the name 
Kolmogorov-Zakharov cascade for wave turbulence). The 
goal of this article is to check directly the hypothesis of 
scale separation in flexion wave turbulence in an elastic 
plate. 

Testing the validity of WTT in experiments is chal- 
lenging because a measurement resolved both in space 
and time is required in order to probe the structure of 
the turbulent wave field. In this respect, the time re- 
solved profilometry technique applied to water surface 
waves and waves on a thin elastic plate is a very valu- 
able tool [1, H|. In the case of water waves it was ob- 
served that the waves develop strongly nonlinear struc- 
tures (crests and bound waves) that explain the disagree- 
ment of the observed wave spectra with the WTT pre- 
dictions [B-Q : for instance the spectral exponents of the 
wave spectrum predicted by the WTT are not observed 
in experiments. The case of waves in a thin elastic plate 
is somewhat different: the space-time spectra of such tur- 
bulence show that the waves are indeed weakly nonlin- 
ear d, IH, H3]: the dispersion relation remains as a sin- 
gle energy line in the (k, uj) space (as opposed to water 
waves) and it is weakly shifted from the linear dispersion 
relation. Despite this agreement with the WTT require- 
ment of weak nonlinearity, the observed wave spectra are 
not in agreement with the theoretical predictions [lll - 
[l3j . This disagreement between experiment and theory 
may be attributed to two effects. First, the forcing is 
localized in space and thus requires some propagation 
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distance to randomize the wave phases and to develop 
the weak turbulence energy cascade This effect im- 
pacts the observed spectrum because the measurement 
region cannot be located far enough from the forcing in 
a region where the stationary energy cascade is effective. 
A second effect may be in action due to the finite size 
of the plate: for finite systems, the modes have discrete 
frequencies. This discreteness may hinder the occurrence 
of resonances if the separation between frequencies is so 
large to prevent resonance conditions to be fulfilled [l4j |. 
Assuming a system of size L, the wavenumbers of the 
neighboring modes should be separated by 5k oc n/L. 
This translates through the group velocity v g in a fre- 
quency separation 5u) sep = v g (uj)5k. Such modes have a 
spectral width Aw due to dissipation (then Aw oc 1/Td 
which is a lower bound of the spectral width) and to 
the nonlinear exchanges of energy that alter the tem- 
poral coherence of the waves. In such a case one ex- 
pects an additional contribution proportional to 1/Tjvl 
which dominates over the dissipation if the scale sepa- 
ration Tjvl < T d is fulfilled. The WTT assumes that 
the frequencies can vary continuously: this is valid if the 
spectral width of the modes is larger than their separa- 
tion (Aw 3> duj se p). The discreteness of the modes can 
be neglected if 1/Tnl 3> Suj sep . If this is not the case, 
the separation of the modes prevents the Kolmogorov- 
Zakharov (KZ) cascade to develop. The dynamics of such 
discrete wave turbulence resembles rather that of a finite 
dimensional nonlinear system [lij ]. In such a system, the 
energy transfer to small scales is strongly affected com- 
pared to the KZ cascade. The continuous KZ cascade is 
expected to be observed in the range of frequencies where 

T < T NL < l/<5w sep . (2) 

Depending of the ordering of Td and 1 / 8ui sep a range of 
discrete turbulence could be observed as well. 

The experimental setup is similar to the one described 
in Q: a stainless plate (lm x 2m x 0.4mm) hangs ver- 
tically, clamped at the top on its short side. An elec- 
tromagnetic shaker is attached slightly above the inferior 
edge of the plate (see fig.^Ih)) driven at 30 Hz with vari- 
ous amplitudes. A videoprojector displays a line pattern 
on the plate. The deformation of the plate yields a defor- 
mation of the pattern. This deformed pattern is recorded 
by an high speed camera at a frame rate up to 10, 000 
frames/s. The deformation field C( r >£] is obtained by a 
2D demodulation of the movies (see [J]). Eventually we 
compute the velocity field v(r, t) = d(/dt. 

Our previous work show that the scale separation T <C 
Td is verified in wave turbulence in an elastic plate [Toj . 
We studied the decay of energy of the wave turbulent 
field after stopping the forcing. A first period of decay is 
associated to the nonlinear cascade and is followed by a 
final exponential decay. We extract from this exponential 
part the dissipative characteristic time Td shown in fig.[T] 
The ratio of Td and the period of the wave is over 100 in 
the accessible range of wavelengths. This scale separation 
enables the existence of an intermediate nonlinear time 
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FIG. 1: (Color online) Squares: Measured dissipative time Td 
versus k/2ir = 1/A. Dashed line: lorentzian fit Td = (0.73 + 
0.025(fc/27r) 2 ) -1 used in the following. Insert: T^w versus 
k/2-K (semilog scale). 

scale. 

In order to check the above scale separation of the 
nonlinear time scale in respect to T and Td, we perform 
an analysis of the dynamics of wave packets through a 
wavelet decomposition. The signal processing yielding 
the estimation of the nonlinear time scale is described in 
part [TTJ The scale separation is discussed in section Mil 

II. DESCRIPTION OF THE WAVELET 
ANALYSIS 

A. Definition of the nonlinear time Tjvl 

In the framework of the WTT, the amplitude of the 
wave is supposed to be slowly modulated by the nonlin- 
ear coupling between the waves. Let us write «k(£) = 

— J v(r,t)e^' r d 2 r the spatial Fourier transform of the 

velocity field. The WTT assumes that v^t) can be writ- 
ten as 

Vk(*) = a k (i)e-^ kt (3) 

where a^(t) is the amplitude which evolves slowly in time 
compared to the period of the wave. We define the non- 
linear time scale as 

where I\(t) = (a,k(t)a^(t + r )) i s the temporal corre- 
lation of the slow modulation of the wave. The average 
( ) is a statistical average on realizations of the turbulent 
field and is estimated as a time average in our case (the 
system is statistically stationnary) . Unfortunately 2?k 
cannot be estimated directly from the experimental data 
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because the velocity field is measured only on a fraction 
of the plate. A wavelet decomposition is performed to 
gain access to the slow modulation correlation function 



at time i on a Gaussian wave packet at position R and 
wavevcctor k. 



C. Correlations of wavelet coefficients 



B. The wavelet family 



a = 15 cm 
A = 5 cm 



k 



FIG. 2: (Color online) Schematics of the plate (figured by the 
vertical rectangular box) of size 2x1 m 2 . Left: Example of 
the Gabor wavelets used in the analysis. The width of the 
Gaussian envelop is a = 15 cm and k/2n = 15 m _1 . Right: 
collection of positions of the center of the wavelets R used in 
our analysis, k is chosen horizontal. The bottom black dot 
corresponds to the position of the electromagnetic shaker. 

In order to extract the dynamics of a wavetrain at some 
wavevector k and position R we use the continuous fam- 
ily of 2D Gabor wavelets: 



Gk,R(r) = 



1 



: exp 



R) 



2cr 2 



exp(jk • r) . (5) 



Such a wavelet is a complex exponential of wavevector k 
modulated by a Gaussian envelop of width a and centered 
in R. The choice of the width of the envelop results from 
a balance between the locality in space and in scale. In 
all the following we use a constant width a = 14.4 cm 
for convenience (rather than fixing the product ak for 
instance). An illustration of the wavelet is shown in fig. [21 
We define the wavelet coefficients as 



tSk,n(t) = J v(r,t)Gl R (r)d 2 



(6) 



where * is the complex conjugation. The wavelet coef- 
ficient is the projection of the turbulent mix of waves 



In order to extract the nonlinear dynamics of the tur- 
bulent field, we investigate the time correlations of the 
wavelet coefficients: 

A,r,r'(t) = (iik,RW^(( + r)) (7) 
Cic,r,r'(t) = ('*2k,R(i)«5*k ) R'(* + T )) (8) 

A is the correlation of the amplitude of a wave packet 
propagating along the same k and between the positions 
R and R' and at times separated by r. C is a similar cor- 
relation but for wave packets propagating along opposite 
wavevectors. Here we study only correlations between 
positions on a line parallel to the horizontal wavevector 
k as shown in fig. [3] The same correlations are computed 
for a plate at rest as well in order to get an estimate of 
the noise contribution. Assuming that the noise and the 
signal are independent, the contribution of the noise can 
be subtracted from the direct estimation of A so that to 
increase the dynamics of the measurement and in par- 
ticular to remove contributions due to the flickering and 
the scanning of the videoprojector (acting at very low k 
and at 60 Hz and its harmonics). 



D. Extraction of the nonlinear dynamics 

In this subsection, we describe how the slow modula- 
tion correlation T> can be computed from the measured 
coefficients A and C. The correlation A can be rewritten 
in Fourier space as 

-Ak,R,R'(r) = j (f; q (t)^,(t+r))G^ R (q)G k , R ,(q')rf 2 qrf 2 q' 

(9) 

where Gk,R(q) is the spatial Fourier transform of the 
wavelet Gk,R(r). 

In the framework of the WTT, two Fourier modes of 
the turbulent field are supposed to be uncorrelated at 
different wavevectors (correlations appear for at least 3 
waves). It has been checked in this experiment for 2- wave 
equal time correlations in [9|. Using this hypothesis, one 
gets: 

(5 q (t)tft(t + r)) - (S q (t)5*(i + T))5(q-q0 (10) 
= £(q-q')2? q (T) e -^ T (11) 

The correlation can thus be rewritten as 
■Ak,R,R'(r)= yG k3 (q)G k , R ,(q)I? q (r) e ^ T d 2 q. (12) 

We assume in the following that 2? q (r) is varying more 
slowly as a function of q than the Fourier transform 
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Gk,R.(q)- By doing this assumption, we suppose that the 
extension of the wavelet (in real space) is narrower than 
the spatial coherence of the wave. As the wavelet has a 
Gaussian envelop, Gk.R(q) is also Gaussian centered on 
q = k. With the above assumption one can neglect the 
variation of T> in the equation (fl"2|) so that: 



^k,R,R'(r)w2J k (T) J G k , R (q)G k , R ,(q)e^ r d 2 q. 

(13) 

The integral in this equation can be calculated using the 
dispersion relation of the waves in the elastic plate to = 
ck 2 and yields 

A ( \ V ^ I (AR-v g (k)r) 2 ^ . 4 
A,R,R'( T ) = i = „ CX P I -■>-■>. I e 



1 



4a 2 (1 



-) 



(14) 

where AR = R' R, v s (k) = 2ck is the group velocity 
and $ = lot + 4>{uj, t) where <f> is a phase factor. 

The correlation A of the wavelet coefficients contains 
the variations of the correlation of the modulation of the 
wave V but altered by the factor 



1 



: exp 



(AR-v ff (k)r)' 



1 



4a 2 (1 



-) 



(15) 



which is due to the dispersion of the wavelet during its 
propagation from R to R' at the group velocity. The 

/ 2 2\~ 1 / 2 

prefactor I 1 + J can be easily compensated. We 
study rather the compensated correlations in the follow- 



ing: A c = yjl + Zg-A and C c = yjl + ^FC. The mod- 
ulus of the compensated coefficient A c is directly linked 
to V: 



IfeWhl^Wlexp 



(AR 



(k)r) 



4^(1 + ^) 



(16) 



Wave packet propagation over a distance AR at the 
group velocity v s account for the Gaussian additionnal 
term. Note that in our configuration, this extra term 
equals 1 for r = e |y A ^)| with e being the sign of AR • v g . 
This time is the flight time of the wave packet at the 
group velocity. Thus by scanning a collection of positions 
R, R' the magnitude of the correlation \T>\ can be recon- 
structed at all time lags r and the nonlinear time scale 
Tnl can be computed. We emphasize that the method 
has been checked to be robust: the result does not de- 
pend on the choice of the wavelet envelope (although we 
restrain our study to Gaussian wavelets in this paper for 
tractability sake) nor on the choice of the width of the 
wavelet a. 



E. Experimental results and rebounds 

In this subsection we present a comparison between 
the experimental data and eq. (I16|) predicted by our 
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FIG. 3: (Color online) Modulus |-4 c (t)| (red squares: ex- 
perimental results; blue solid lines: analytical predictions) 
and real part jy4 c (t)| (red circles: experimental results; blue 
dashed lines: analytical predictions) of the compensated cor- 
relation functions -4 k R R+AR (f). N° adjustable parameter 
is used, (a) k/{2ir) = 10.6 m" 1 , AR = 0; (b) k/(2n) = 
20.2 m~\ AR = 0; (c) fc/(27r) = 20.2 m" 1 , AR = 0.27 m. 



model: examples of compensated correlation functions 
r R' ( T ) an d C£ R R , (r) are shown in fig [3] and fig. [4] 
Figure [3] compares measured correlations for short times 
with the expression given in eq (|16j) . The modulus of 
the correlations at the same positions R = R' (fig[3Ja) 
and (b)) exhibit a maximum for r = 0. As r grows, the 
wavelet propagates away from its original position and 
this modulus decreases to in the Gaussian fashion pre- 
dicted by equation (Ti"4"l) . The real part oscillates as pre- 
dicted by the term exp (J [cot + <fi]). As expected, a wave 
at a larger wavenumber propagates further and oscillates 
faster than for smaller wavevectors. An example of cor- 
relation for different positions is displayed on figure OD^c) . 
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FIG. 4: (Color online) Example of correlations at k = 
16.0 m~ . The real part is shown as the oscillating thin red 
line and the modulus as the thick black line, (a) ■Ak,R,R.( T ) 
for R at the center of the plate. The vertical dashed lines 



show the times 2nL/v g for n 



-2 to 2. (b) Ai 



(r) 



for AR = 0.27 m. The vertical dashed lines show the times 
(2nL + AR)/v g for n = -2 to 2. (c) C£ jB . iR (r) for R at 
the center of the plate. The vertical dashed lines show the 
times (2n + l)L/v g for n = —3 to 2. In the three cases, the 
correlations have been normalized by -4k,R,R,(0). 



The maximum is delayed due to the flight time. Mea- 
surements shows a very good agreement with predictions, 
although a small departure from the expected Gaussian 
envelop is observed for small wavevectors / large wave- 
lengths (figure [3ja)). Indeed, the second statement of 
our derivation is not met for small wavenumbers: the 
wavelength and the size of the wavelets are comparable. 

The measured correlation coefficients are plotted for 
larger time lags on figure |3J The autocorrelation coef- 
ficients A (fig HJa,b)) exhibit some secondary maxima. 
They correspond to a wave packet that has bounced on 
the borders of the plate. After one rebound, the wave 
packet propagates back with the opposite wavevector — k. 
Thus no correlation is observed until it bounces a second 



time and propagates again with its original wavevector. 
Correlation reappear as it comes back to its initial po- 
sition R. These reflexions account for repeated maxima 
with a delay Tfl y = 2L/v g corresponding to the wave 
packet round trip accross the plate. These rebounds can 
be used to extend the range of accessible values of the 
time lag for the reconstruction of T>. The magnitude of 
the reflexion coefficient at the border is expected to be 
very close to unity for a free border moving in air (whose 
density is about 4 orders of magnitude lower than steel) . 
The relative amplitude of the maxima is decaying with 
the time lag because of the decoherence of the wave as it 
propagates. 

Figure (DJc) shows C c for R = R'. As expected, the 
correlation is zero at small r because the wave packet 
has not enough time to propagate to the border of the 
plate and back to the orginal position. A maximum of 
the modulus of the correlation is observed for a time lag 
t = L/vg. This corresponds to the time needed for a 
wave packet to travel to the border, bounce on it and 
come back to its original position counterpropagating. 
Secondary maxima are also observed with a period Tfl y . 
By collecting all values of all maxima of A c and C c one 
can build the modulus of T> for quite long time lags (until 
the signal falls under the noise level). 
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FIG. 5: (Color online) Modulus of the correlations at the same 
position R = R' as a function of time and wavevector. Colors 
are log coded, (a) Autocorrelations |.4 C |. Green (light gray) 
dashed lines: flight times for 2n rebounds t^ y — 2nL/(v g (k)); 
black dashed line: flight time over the width of the wave 
packet (b) Crosscorrelations \C C \. Green (light gray) dashed 
line: flight times for 2n— 1 rebounds t„ y = (2n— l)L/(v g (k)). 
At each value of k, the correlations have been normalized by 
At 

r.r(0) so that the modulus of the correlation is between 
and 1. 



Figure [5] presents in its upper part the correlation 
R r (t) as a function of both k and r and C£ R r (t) 
in the lower part. The dashed lines correspond to flight 
times over 2nL for A or (2n + 1)L for C at the group 
velocity. The local maxima of correlation follow these 
lines: this validates our interpretation of the secondary 
maxima as due to the rebounds. 
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FIG. 6: (Color online) (a) Color plot of 

j^k r,r'(*)| / |"4k,R,R,(0)| versus time t and position 
R'. Wavcnumbcr k/2iv = 13 m _1 ; the position R is at 
center of the plate. Colors are log-coded, (b) Crosscorre- 
lation |C£ R R / (t) j / |-4k,R,R(0) | (c) Combinaison of the two 

functions A c and C c : T = J \A C \ 2 + |C C | 2 . 



A picture of the travel of the wave packet can be build 
by plotting the modulus of the correlation in a space time 
representation as in fig. [6l Figure EJa) and (b) show the 
modulus of A c and C c respectively. Figure displays 
only the wave packet traveling from left to right whereas 
fig. (njb) shows only the reverse motion. The trajectories 
appear slightly deformed in a S shape due to the fact 
that close to the border of the plate, the wave packet is 
partially reflected in the other direction and the wavelet 
is truncated by the border. The two functions can be 
combined to show the full trajectory in fig. EJc). The 
correlation is seen to decay in magnitude as time goes by 
along the trajectory corresponding to the loss of temporal 
coherence of the wave packet due to the energy exchanges 
with the other waves by the nonlinear coupling. 

Figure [7] shows the space-time representation of the 
decay of the correlation along the trajectory of the wave 
packet for 4 different values of the forcing. As the forc- 
ing is increased, the number of rebounds that can be 
observed is also decaying: 6 bounces can be seen in the 
60 ms time window at the lowest forcing while only three 
are visible at the strongest forcing intensity. The stronger 
the forcing, the higher the nonlinearity, the faster the en- 
ergy transfer and thus the smaller the nonlinear time. 
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Position [m] Position [m] 

FIG. 7: (Color online) Color plot of T (see previous figure) 
versus time t and positions R'. Wavenumber k/2-K — 13 m _1 ; 
the position R is the center of the plate. Colors are log-coded. 
The forcing power is P , 7P , 27P and 59P (from (a) to (d)) 



F. Nonlinear time measurement 
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FIG. 8: (Color online) Superimposition of the modulii 

\AZ 

rr'( t )I (normalized by |-4kR R (0)| at the center) at a 
given k and for a collection of positions R' (excluding posi- 
tions too close to the border for which the wavelet can be 
truncated) . R is chosen at the center of the plate. The thick 
line corresponds to R = R'. Black thin lines: -4kRR'( T )i 
cyan (light gray) thin lines: C£ RH /(t). Insert: squares cor- 
respond to the value of the local maximum. The line is an 
exponential fit of the decay. 

In order to rebuild the modulus of the correlation 
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2?k(r), one should collect the values of the local max- 
ima of the correlations A c and C c for various pairs of 
positions R, R'. An example of collection of all such cor- 
relation curves is shown in fig. [5] for given k and forcing 
intensity. The maxima are extracted by a local Gaussian 
fit and shown in the insert. The decay of the maxima 
with the time lag is seen to be exponential so that 



|Pk(r)| 
l^(0)| 



cx exp - 



T NL (k) 



(17) 



which permits a direct extraction of the nonlinear time 
scale. 



III. DISCUSSION 



A. The nonlinear time scale 




FIG. 9: (Color online) Evolution of the nonlinear time scale as 
a function of the wavenumber for different forcing. Symbols 
are experimental data for different values of the forcing inten- 
sity (intensity increasing downward): Po, 7Po, 27Po, 59Po, 
lOOPo, 149Po, and 207Po. This colorchart for injected powers 
is used through this section for other figures. 

The evolution of the nonlinear time scale with the 
wavenumber and the forcing intensity is shown in fig. 
The estimation of Tnl was not performed for fc/2-7r less 
than 5 as the wavelength is too large to fit in the Gaussian 
modulation of the wavelet (the wavelength is compara- 
ble to a). At large k, the estimation is performed until 
the signal over noise ratio gets too degraded to allow 
any computation of the correlations. Thus the highest 
wavenumber that is reached depends on the amount of 
data gathered at the considered value of the forcing as 
well as on the frame rate of the camera (the highest the 
frame rate, the lowest the amount of light, the highest 
the noise level). The nonlinear time may be considered 
as a measure of nonlinearities: in this respect, we expect 
this time to diverge for linear problems and to decrease 



as the nonlinearities grows. Figure [9] shows that Tml is 
indeed a growing function of k and a decreasing function 
of the forcing P. The evolution of T^l with k is affine. 
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FIG. 10: (Color online) Rescaled nonlinear time TnlP ' 5 as 
a function of the rescaled wavevector fcP _1//6 for different in- 
jected power P. 



A cutoff frequency was previously measured from the 
wave spectrum 13]: /* oc P 1 / 3 which translates through 
the dispersion relation into k* oc P 1 / 6 . This cutoff fre- 
quency corresponds to the transition to a dissipativc 
regime. Here we look for the following rescaling of the 
nonlinear time scale: 



T NL (k)acP- a f(kP-V* 



(18) 



The curves for different forcing collapse on a master curve 
when using a — 0.5 (fig. ITD|) . This scaling is similar to 
that previously observed for the power spectrum [l2], [H| ■ 
In the framework of the WTT, one would expect that 
Tnl scales as P -1 / 3 for 4-wave resonance as predicted by 
the theory [ll| or as P" 1 / 2 for 3- wave resonance. Here, 
the rescaled time is affine with the rescaled wavenumber 
so that a pure scaling is not observed. The nonlinear 
time follows a scaling law 



T NL = /3P~ 1/2 + jkp- 2/3 



(19) 



with constant /3, 7. 



B. Time scale separation 

The experimental time scales are compared with the 
dissipative time scales and the period of the wave in 
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FIG. 11: (Color online) Comparison of the nonlinear time 
scale with other relevant time scales. Symbols are experimen- 
tal data for Tnl at different values of the forcing intensity. 
Upper dashed line: twice the dissipative time Td measured by 
the decay of energy in the unforced case (see fig[TJ. Middle 
dashed line: 2/Suu 3ep = L/(irkc). Lower dashed line: 1/w. 
The forcing is operated at 30 Hz and is shown with the verti- 
cal dash-dotted line. 



fig. 111! The figure displays rather twice the time Td es- 
timated from the decay of the energy in the decline case 
as the energy is quadratic in the amplitude of the wave 
thus the decay time of the amplitude is twice that of 
the energy. The figure shows 1/oj rather than the pe- 
riod 2it/uj for the following reason: the correlation of the 
magnitude of the wave is observed to be exponentially de- 
caying. Its Fourier transform follows thus a Lorentzian 

shape oc 1/ ( uj 2 + -fk— ) • Thus Tn l should be compared 

V 1 nlJ 

directly to 1/uj. 

The product Tn luj is shown in figure I12f a) . The non- 
linear time scale is seen to be over an order of magnitude 
above w for most of the accessible wavenumber range so 
that the scale separation loTnl 3> 1 is verified. At the 
highest forcing intensity the ratio is only a factor 2 at 
the forcing scale suggesting that the hypothesis of weak 
nonlinearity may be disputable at this scale. 

The ratio 2T d /T NL is displayed in figure [T^b) . The 
nonlinear time scale is over two orders of magnitude 
shorter than the dissipative time at the forcing scale. The 
ratio is always larger than 3 at the highest wavenumber 
that can be analyzed (which is always above the observed 
dissipative cutoff). In the framework of the Weak Turbu- 
lence Theory one would expect the conservative energy 
cascade to operate until Tnl — Td and to be followed by 
a dissipative regime. Here it appears that the systems 
reaches the dissipative regime at wavenumbers for which 
the nonlinear time are still larger than the dissipative 
times. The cascade is stopped by another phenomenon 
before reaching the expected dissipative region. This ob- 
servation is attributed to finite size effects in the following 
section. 
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FIG. 12: (Color online) Ratios between different characteris- 
tic times as a function of k for various values of the injected 
power: (a) loT nl ; (b) 2T d /T NL ; (c) T NL SLu Bep /2 



C. Finite size effects and transition to frozen 
turbulence 

Finite size effects can be neglected when the half sep- 
aration of the modes 8uj sep /2 is larger than 1/Tnl- The 
quantity 2/5u) sep is compared with Tnl in fig- 1111 and the 
ratio is shown in figure[TJJc). At low forcing intensity, the 
product TnlSloscp is always larger than one so that finite 
size effects are expected to be dominant over the whole 
interval of accessible wavenumbers. When increasing the 
forcing, the product gets smaller. At the strongest forc- 
ing the ratio is lower than one over the whole interval of 
k so that a behavior close to WTT should be expected. 
At intermediate values of the forcing a crossover is seen 
from a regime were finite size effects should be negligible 
(at low k) to a regime dominated by finite size effects (at 
large k). The crossover region corresponds to laminated 
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turbulence as described in EMI. In this regime, the 
resonance conditions on frequencies uji + w 2 = w 3 + w 4 
(expected for 4- wave resonance [TT}| ) are made less nu- 
merous due to the discretization of the modes. Nonlinear 
spectral widening allows for approximate resonances if 
the nonlinearities are strong enough so that a continuous 
regime can be observed. In the opposite case, only a few 
clusters of frequencies can be resonant and the cascade 
operates with a reduced efficiency. If the nonlinearity is 
too weak, the cascade is frozen when Ti 



NL 



Slj s 



so that 

the energy cannot be transferred to higher wavenumbers 
and is ultimately dissipated. In this case, the energy 
cascade stops before reaching the dissipative regime ex- 
pected from WTT (at which T NL ~T d ). 



U 



^ 10 1 



|4 1 



J 



(a) 



i t-** 



\l m ***** * 



5 10 15 20 25 30 

k/(2n) [m" 1 ] 




1000 1500 2000 

/ [Hz] 



FIG. 13: (Color online) (a) Symbols: measured Tnl- Dashed 
line: 2/Susep- Red filled squares: crossover wavenumber k c 
given by TVlCO = 2/5u sep (k c ). (b) Solid lines: Stationary 
power spectrum of the velocity extracted along the dispersion 
relation and for k x = (see text). Forcing is increased from 
bottom to top. Plain squares: crossover frequency f c reported 
from (a), (c) Synthetic spectrum incorporating the nonlinear 
widening of the modes (see text). The width of the modes is 
1/(2ttTnl) interpolated as a function of frequency from (a). 
Curves are vertically shifted for clarity on fig (b,c). 



Figure [13] analyses the experimental observation of the 
widening of the modes due to nonlinear effects. The 
cross-over wavevector k c defined as: 



T N h{k c ) = 



Suj sep (k c ) 



(20) 



is first determined on fig. [i~37a) for different injected 
power. Using the dispersion relation, this crossover 
wavenumber is associated with a crossover frequency / c : 



f c = Ck 2 c /(27T) 



(21) 



By performing Fourier transforms both in time 
and space, one computes the power spectrum density 
E(k,u>) = (|w(k, w)| 2 ). Figure \W(b) shows a part of the 
power spectrum density E{k x = 0,k y = ^/uj/c,lu), i.e. 
along the dispersion relation in a cut k x — 0. A tran- 
sition can be observed between 2 distinct regimes: the 
spectrum is smooth and continuous at low frequencies 
whereas it displays some clearly separated peaks at high 
frequencies. This transition is well described by our cri- 
teria (eq (|2ip ): the crossover frequency extracted from 
figure fT3Ta) separates the two distinct regimes as shown 
in figure fl~3T b). A toy model can be build to incorporate 
the nonlinear widening of discrete modes. The frequency 
of the modes are given by /„ = c(mr / L) 2 / (2tt) and each 
mode is broadened following a Lorentzian variation of 
width 1/(27tTjvl). A density of energy can be build in 
the following way: 



p(/)«E 



i 



[2nT NL (f - f n )Y + 1 



(22) 



This synthetic density of mode p(f) is shown in fig.[T37c). 
The overall characteristics of fig-ITBTb) are recovered (ex- 
cept for the general decay which is not taken into account 
in the model). 



D. Width of the dispersion relation 

Our previous work demonstrated that the energy was 
indeed localized in the space-time power spectrum den- 
sity ECk,ui) in the vicinity of the linear dispersion rela- 
tion [J, [9J. Due to weak interaction between waves, the 
observed dispersion relation is slightly shifted from the 
linear one. The observed dispersion relation shows also 
a finite width Aw that can be attributed to two distinct 
causes: the nonlinear widening Aujnl = I/TWl as ex- 
pected from the WTT and the finite resolution due to the 
finite size of the measurement picture Acj W i n = 2cfcAfc = 
2ckir/l. These phenomenon are independent. Hence the 
effective width is expected to be close to: 



(Aw) 2 ~ (Awtvl) 2 + (Aw win ) 2 



A 

1 NL 



+ Bk 2 



(23) 



where A and B have been substituted to the above the- 
oretical value to account for geometrical constants in 
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FIG. 14: (Color online) (a) Circles: (Auj/k) 2 as a function 
of (Tjvifc)" 2 for different injected powers. Dashed line (eye- 
guide): 0.6a: + 25. (b)A(T NL k)~ 2 + B (circles) and (Aco/k) 2 
(squares) plotted as a function of k for one given injected 
power P = 149P - A = 0.6, B = 25. Solid line: value of B. 



the experimental evaluation of Aw and in the effective 
shape of the energy profiles due to both phenomena. 
A is expected to be close to 1 and B of the order of 
(2irc/l) 2 ~ 65. The integrated over angles power spec- 
trum E(k,u) is fitted with a local Gaussian along u> at 
given k to measure the width Aw(fc) (as in [9(). We 
compare this width with the nonlinear time measured 
from the wavelet analysis presented in this paper. Ac- 
cording to equation (|23[) . a linear relation is expected be- 
tween (Aco/k) 2 and (Tjvxfc) -2 - This is indeed observed 
in fig. [T47a). The extracted values A and B are used to 
compare Aw and Tnl as a function of k in fig. 114( b) and 



the agreement is very good. This observation validates 
the fact that the nonlinear time is responsible for the 
widening of the dispersion relation as predicted by the 
WTT (even if the observed scaling of Tml is not the one 
predicted by the WTT). 

IV. CONCLUSION 



We proceeded to a direct extraction of the nonlinear 
time scale at various wavenumbers and forcing intensities 
by studying the temporal correlations of wavelet coeffi- 
cients. With this technique, we have checked that the 
scale separation hypothesis Td 3> T^l 3> T is reason- 
ably valid in our system which is indeed weakly nonlin- 
ear. This hypothesis is one the main requirements for 
the validity of the Weak Turbulence Theory. Neverthe- 
less, the scaling of the power spectrum were previously 
observed to be in disagreement with the theory. It ap- 
pears also that the cascade does not really proceed to the 
expected dissipative regime for which Tn l is comparable 
to the dissipative time scale. This is due to the fact that 
the second major hypothesis of the theory is not valid: 
the asymptotically large system. Finite size effects lead 
to discrete modes. This discreteness can be ignored if the 
spectral widening due to nonlinear energy transfers over- 
comes the frequency separation of adjacent modes in the 
k space. We checked that this requirement is not fulfilled 
at weak forcing whereas it is valid at low wavenumber for 
strong forcing intensities. One would expect a range of 
large wavelength for which the WTT should be valid. 
But this range is most likely altered by another sort of 
finite size effects related to the forcing and previously 
analyzed [Icl |: in the vicinity of the forcing point, the 
phases of the waves with large wavelengths need some 
distance to get randomized. The plate is most likely too 
small for the measurement region to be far enough from 
the forcing so that to really observe a pure WTT regime. 
Nevertheless, many qualitative feature of the Weak Tur- 
bulence Theory and of laminated turbulence have been 
observed in this system due to the possibility to perform 
advanced space-time analyses of the turbulent field made 
possible by the time resolved profilometry technique p|. 
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